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Abstract 

P^ I We generalize the results of Ambrosio [Invent. Math. 158 (2004), 227-260] on the exis- 

Ph ' tence, uniqueness and stability of regular Lagrangian flows of ordinary differential equations 

(-H I to Stratonovich stochastic differential equations with BV drift coefficients. Then we con- 

struct an explicit solution to the corresponding stochastic transport equation in terms of the 
stochastic flow. The approximate differentiability of the flow is also studied when the drift 
is a Sobolcv vector field. 
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^ 1 Introduction 

O 

Q I Let Aq, Ai,- ■ ■ , Am be vector fields on M and wt = {wj:, • • • , w^^) an yn-dimensional standard 

Brownian motion defined on a probability space (fi, J-", P). Consider the Stratonovich stochastic 

differential equation (abbreviated as SDE) 



dXt = J2 MXt) o dwi + AoiXt) dt, Xo = X. (1.1) 

1=1 

It is well known that if At G C^2+'^(R'^,R^) (f = I,--- ,m) and Aq G CjJ+''(R^,R^), then the 
above equation has a unique solution which defines a stochastic flow of C^-diffeomorphisms on 
M . Here C^"*" means that the vector field ^o and its first order derivatives are bounded, and 
the derivatives are Holder continuous of order 5 > 0. These conditions on the boundedness 
of the vector fields and their derivatives were relaxed in [18], by allowing the local Lipschitz 
constants on the balls centered at the origin to grow as fast as the logarithmic function. In 
the case (5 = 0, it is proved in [12] that under the same growth conditions, (1.1) still gives 
rise to a flow of homeomorphisms on R*^. This result is generalized in [13] to the case where 
the drift coefficient ^o satisfies only the general Osgood condition, at the price of the C^, 
regularity for the diffusion coefficients; if in addition the distributional divergence of ^o exists 
and is bounded, then the Lebesgue measure is quasi-invariant under the action of the stochastic 
flow of homeomorphisms (cf. [19]). 
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On the other hand, the ordinary differential equation (abbreviated as ODE) 

dXt = Ao{Xt)dt, Xo = x (1.2) 

with Sobolev or even BV coefficient has been studied intensively in the last three decades. Here 
^0 can be a time-dependent vector field. The existence of quasi-invariant flow of measurable 
maps associated to a vector field Aq with Sobolev regularity was first studied by Cruzeiro [8]. 
A thorough treatment was carried out by DiPerna and Lions in the celebrated paper [9], where 
the authors deduced the existence and uniqueness of fiows generated by (1.2) from the well 
posedness of the corresponding transport equations. Similar results were obtained in [6] by 
taking the standard Gaussian measure as the reference measure. Ambrosio [1] generalized the 
results to the case where ^o has only BV regularity by considering the continuity equation. S. 
Fang [11] gave a short introduction to the theories mentioned above. The extension of these 
results to the infinite dimensional Wiener space have been done in [3, 14]. Using the local 
maximal function, Crippa and De Lellis obtained in [7] some new estimates which allow them 
to give a direct proof of the existence and uniqueness of the DiPerna-Lions fiow. 

Inspired by these studies of ODE, there have been several attempts to solve the SDE with 
Sobolev coefficients. Following the method in [7], X. Zhang [23] showed the existence and unique- 
ness of the stochastic flow of measurable maps generated by Ito SDE with Sobolev coefficients, 
provided that the derivatives of the diffusion vector fields are bounded. The SDE with BV drift 
vector field was also considered in this paper, but the diffusion coefficients were assumed to be 
constant. In [15], the authors took the standard Gaussian measure as the reference measure 
and proved a priori estimate on the L^ norm of the density of the flow, which enabled them to 
construct the unique flow associated to (1.1), provided that the gradients of the diffusion coef- 
ficients and the divergences with respect to the Gaussian measure are exponentially integrable. 
In the recent work [24], X. Zhang studied the Stratonovich SDE with drift coefficient belonging 
to W^q'c (M , M ), and he also provided a Freidlin-Wentzell type large deviation estimate for the 
stochastic fiow. 

In the present paper we consider the Stratonovich SDE (1.1) with BV drift vector field. 
Our method is based on Ocone-Pardoux's decomposition [21] of the fiow generated by (1.1) into 
the stochastic fiow of the diffusion part, and a flow associated to random ODEs whose driving 
vector field is a transformation of the drift coefficient Aq by the stochastic fiow. This approach 
was applied in [13] to deal with the Stratonovich SDE with drift satisfying the general Osgood 
condition. Our main result can be stated as follows {Cd is the Lebesgue measure on 
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Theorem 1.1. Assume that Ai,--- ,Am G C^~^''(M'^,M'^), and the drift Aq satisfies 

(1) Aq has sublinear growth, i.e. \Aq{x)\ < C(l + |x|^~^"), x € M ; 

(2) Ao€BVioe(K^]R'^); 

(3) the divergence D ■ Aq <^ Cd with locally hounded density function d\Y{AQ). 

Then the equation (1.1) generates a unique stochastic flow of measurable maps Xj : M — t- M 
which leaves the Lebesgue measure Cd quasi-invariant. 

This result will be proved in Section 3. Note that the sublinear growth of Aq ensures that 
the vector field ^o defined in (3.1) has similar growth (Lemma 3.1), which in turn implies the 
classical growth estimates on the solution of the ODE. The sublinear growth of ^40 also allows 
us to assume only the local boundedness of the divergence div(^o) (as in (ii) of Theorem 6.2 
in [1]), compared to the global boundedness required in [2, 9]. The starting point of the proof 
of Theorem 1.1 is the relationship between the distributional derivative (a Radon measure in 



the present case) of the drift Aq and that of the transformed vector field ^o defined in (3.1). 
This will be done in Lemma 2.2, where we also show that if the divergence D ■ Aq is absolutely 
continuous with respect to Cd, then so is Z? • ^Oj and they are related to each other by a quite 
simple equality. 

We study in Section 4 the stability of the solutions to (1.1) when a sequence of vector fields 
Aq converge in some sense to Aq. For this purpose, denote by B{R) the ball centered at the 
origin with radius R. 

Theorem 1.2. Assume the conditions of Theorem 1.1. Let Aq : W^ ^ W^ he vector fields 
satisfying 

(1) there is C > and Eq G (0, 1), such that sup„>i |^o(x)| < C(l + |x|i-^o), x G M'^; 



(2) Aq converges to Aq in Lj^ 
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(3) for any n > 1, V^q is locally bounded; 

(4) for any R>0, sup„>i ||div(^^)||ioo(5(^)) < +cx). 

Let X^ be the flow associated to (1.1) with Aq being replaced by Aq. Then for any p > 1 and 
T, R> 0, the following convergence holds almost surely and in LP{Q,F): 



[ sup |Xj"(x) - Xt{x)\dx = 0. 

Jb(R) tGfO.Tl 



lim 

It is well known that when the coefficients are smooth, the solution to the corresponding 
stochastic transport equation can be explicitly expressed in terms of the flow generated by 
(1.1), see [13] Theorem 5.1 for the case where ^o satisfies the general Osgood condition. As an 
application of the above stability result, we will show that similar representation still holds even 
in the situation of BV drift. 

Finally in Section 5 we consider slightly more regular drift coefficient Aq G Wj^^^ , showing 
that almost surely, the stochastic flow Xt is approximately differentiable on M*^. This generalizes 
the results in [4, 7] to the stochastic context. 

2 Preparations and known results of ODE 

In this section we give some preliminary results needed in the subsequent sections. Here is the 
definition of BV functions. 

Definition 2.1. A locally integrable function 6 : M'^ — ?> M'" is said to be of class BV/oc tf there 



IS a 
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-valued Radon measure Db = (i?ifr')i<i<d, i<j<-m such that 

^d{D,b>) = - I Vdiijdx, for all ^ G Cli^!^). 

If 6 G C^, then we still denote by Db the function V6 on M . li m = d, we denote by 
D ■ b = tr(D6) = X]j=i-^«^* ^^^ divergence of the BV vector field b, which is again a Radon 
measure on M . In the following det(-) means the determinant of a matrix. For a measurable 
map if :W^ ^ W^ and a Radon measure fi on W^, ip^fi denotes the push forward of the measure 
^hy if (or the "distribution" of ip under fi). Now we prove 

Lemma 2.2. Let b : M.'^ —^ W^ be a BVjoc vector field and if : W^ —?■ M.'^ a C"^ -diffeomorphism, 
then 



(1) the composition ho tp is still a BV^oc vector field and 

D{ho^) = \det{J^)\-^j;[{^-^)#Dh], 

where J^ is the Jacobi matrix of if and J* is its transpose, J* [(c/?^"'^)^Z?6] is the matrix 
product of J* and {ip~^)^Dh; 

(2) if the divergence D -b <^ Cd with density function div(6), then 

D ■[J-\bo^)] = {d\v{J~^),ho^) +dw{b)o^, 
where div(J~-'^) is a vector field whose components are the divergences of the column vectors 

of J:,'. 

Proof. (1) For every i € {1, • • • ,d}, we only have to show that 6* o y? is a BV function. By 
Theorems 2 and 3 in Section 5.2 of [10], there exists a sequence of functions 6^ G BV;ocnC°°(M'^) 
such that 6^ — )■ 6* in Lj^^{M.'^) and the vector valued measures DK^ converges weakly to Db^ as 
n — >• oo. Note that the composition b\^o Lp belongs to C^, hence for any ijj G C^iW^), by the 
integration by parts formula, 



V' [D{1^^ o^)]dx = - I {bi o (^) W dx. 



(2.1) 



We have by the chain rule, -0(6^ o ip) = J* [(D6^) o ip\. It follows from the formula of changing 
variables that 



ij[D{b\o^)]dx= ijr^[{Db\)o^]dx 

[(V^j;)o(^-i)](Z)6^)|det(J^-i)|da; 
I 

as n ^ oo, due to the weak convergence of (Dbl^) dx to Db^. Now we have 

[{^PJ;) o ((^-1)] . |det(J^-i)| d(Z)6^) = / ^ V|det(J^-0 o ^\j; d[iip-%Db- 

^P\det{J^)r'J;d[{^~%Db^], 



J 



where the last equality follows from J^-i o p = (J^) ^. Therefore 

hm / ^lJ[D{bl,o^)]dx= f i;\detiJ^)\-'j;d[{^-%DV]. (2.2) 

Now we consider the limit of the right hand side of (2.1). Again by changing variables, 
- / (6>(^)VV'dx = - / 6J,.[(V^)o(^-i].|det(Vi)|dx 
which converges to 

b' ■ [(V^) o (^-1] . |det( J -1 ) j dx = - / (6* o (^) W dx. 



This combines with (2.1) and (2.2) leads to 

V|det(J^)rij;d[(v9-i)#m*] = - / (6'o^)V^dx, 

which means that 6* o (^ is a BV function and 

D{b^o^) = \det{J^)r^j;[{^-%Db']. 

(2) By the chain rule and (1), it is easy to know that J^^{b o (^) is a BV/oc vector field, and 

D ■ [J-\b o V9)] = (div(J-i), 6 o V9) + {J-\D{b o ^)), (2.3) 

where the second (•, •) is the inner product of matrices regarded as elements in M . By the 
expression in (1), we have 

{J-\D{bo^)) = |det(J^)|-itr[((/,-i)#Z?6] = |det(J^)|-i((^-i)#tr(D6). (2.4) 

Since tr(L>6) = D-b = diY{b)Cd, for any V G C^{W^), 

i)d[{ip-^)#iT{Db)]= [ iij o ^~^) dMDb)] = [ (V^ o c^-i) . div(6) dx 

Tp- [d\Y{b)oip\ ■ |det(J<^)|dx. 

Therefore {ip~^)^iT{Db) = [div(6) o c^] • |det(J<^)|. Combining this equality with (2.3) and (2.4), 
we complete the proof. D 

The next technical result will be used in Section 4. 

Lemma 2.3. Let c^ : R"' — )• M"' 6e a C"^ -dijjeomorphism and p := '^^ '' the Radon- Nikodym 
density function. Then we have 

(p-iVp)o^ = div(j;i). 

Proof. It is well known that p o ip = \det(^J^^^\. Since ip is a C^-diffeomorphism of M , 
the function M"^ 9 x i-> det(J~^)(x) does not change sign. Without loss of generality, we may 
assume that det(j~^) > on the whole M'^, thus pop = det(j~-'^). As V{po(p) = J*[(V/j) o p], 
we have 

(Vp) o p = {J~^rV{p o p) = ( J-i)*Vdet( J-i) . 

Note that Vdet(J~-^) = — [det(J,^)]~^Vdet(J^), so the equality that we should prove is 

Vdet(J^) = -det(J^)j;div(j-^). 

For simplification of the notations we write J = J^ and K = Jz^, then Jij = djp'^, 1 < i,j < d. 
We consider the determinant det(-) as a smooth function defined on R . By the chain rule 
and Jacobi's formula (see [20] Part Three, Section 8.3), we have for any / G {1, • • • ,d}, 

didetiJ) = Y, T^iJ) ■ 9iJij = E det(J)ir,i%<^\ (2.5) 



For any I < j < d, it holds 5ji = Yl'i=i KjiJu = Y.'i=i Kjidiif\ therefore 

d d 

that is, J2i=i ^jii^jiV^) = — J2i=i{9j I^ji)idi^^) ■ Combining this with (2.5) leads to 

d d d d 

didei{J) = -det(J) Y, J2^djK,i){div') = -det(J) ^idi^') Y,{djK,i) 

j=l i=l i=l j=l 

d 

= -det(J) Y, JiidHK-i) = -det(J)(J*div(K))z. 

The proof is complete. D 

Now we recall the definition of the regular Lagrangian flow associated to a time-dependent 
vector field bt (see [2, 7]). 

Definition 2.4. Let b € Ll^^{[0,T] x W^,W^). We call a map Y : [0,T] x R"^ ^ W^ a regular 
Lagrangian flow for the vector field b if 

(1) for a.e. x G M'^, the map [0, T] 3 t ^- Yt{x) is an absolutely continuous integral solution of 

dYt = bt{Yt)dt, Yo = x; (2.6) 

(2) (Yt)#Cd<:/:dforalltG[0,T]. 

Note that this definition is slightly different from that in [7]: in condition (2) we do not 
require that {Yt)#Cd is dominated by CCd on the whole M'^, where C > is a constant. 

Given a measurable map y : R"^ — )■ M'^, we say that Z : R*^ — )• M"^ is a measurable inverse 
map of y if Z is measurable and for a.e. x G R , x = Y{Z{x)) = Z{Y{x)). We often denote by 
Y^^ the measurable inverse map of Y (see [23] Lemma 3.4 for a characterization of this notion). 
In the following theorem we summarize the results concerning the existence and uniqueness of 
regular Lagrangian flow generated by a BV vector field. 

Theorem 2.5. Let 6^ : R — > R be a time- dependent vector field satisfying: 

(1) M£^GL-([0,r]xR'^); 

(2) bt e BVzoc(IR'') for a.e. t € [0,r], and for any R>0, \Dbt\{B{R)) € Ll^{0,T) and 

T 

||div(6t)||ioo(B(^))dt < +00. 

Then the vector field b generates a unique regular Lagrangian flow {Yt : < t < T} . Moreover 
for any t £ [0,T], Yj has a measurable inverse map Y^T : R — ^ R and (Y^~ )#^d = Pt^d with 

pt{x) = exp f / div{bs){Ys{x)) ds j . 

Proof. The first part of this theorem was first proved in [1] for bounded vector field b, and 
then in [2] for the general case (see the remark at the end of Section 5), while the second part 
was proved in [23] for a BV vector field b independent of time (just let the diffusion coefficients 
a be in Theorem 2.6), but the proof for the time-dependent case is similar. D 



Remark 2.6. For the density function pi := d/f '' i ^^ have 

pt{x) = exp(-f dw{bs)[Ys{Yf\x))] ds' 
See [6] Theorem 2.1 where the expressions are given using double time parameters. 

3 Existence and uniqueness of (11) with BV drift 

In this section we prove Theorem 1.1. First we introduce Ocone and Pardoux's decomposition 
of the Stratonovich SDK (1.1) (see [21] PART II or [13] Section 2). Consider the foUowing 
Stratonovich SDE without drift: 

m 

dXt = Y,MXt)odwi, Xo = x. 

i=l 

It is well-known that under the conditions that Ai,- ■ ■ , Am G C'^"'' for some 6 > 0, the solutions 
of the above SDE admit a version Xt{x,w) such that there exists a full subset Oq, for each 
w; G Oq and each t > 0, x — >■ Xt{x,w) is a C^-diffeomorphism of W^. Set ipt{x) = Xt{x,w). Let 
Jt{x) = dx(pt{x) be the Jacobian matrix of ipt : M."^ ^ M'^ and Kt{x) = (Jf(x))~^ the inverse of 
Jt{x). Define for u; € JIq, 

AQ{t,x) = Kt{x)Ao{ipt{x)). (3.1) 

We consider the differential equation 

m = Ao{t,Yt)dt, Y^ = x. (3.2) 

Then the solutions of (1.1) can be expressed as (at least when ^o is smooth) 

Xt{x) = n{yt{x)). (3.3) 

Therefore it is sufficient to study the well posedness of the random ODE (3.2) under the as- 
sumptions on Aq in Theorem 1.1, and then show that the representation (3.3) indeed gives the 
flow associated to the original Stratonovich SDE (1.1). 

In the next lemma, we collect various growth results concerning the stochastic flow ^pt and 
its derivatives, and the random vector field Aq defined in (3.1) (see [13] Lemma 2.2 for a proof). 

Lemma 3.1. Assume the conditions of Theorem 1.1, we have for any T > 0, 

(1) for any a > 1 and P > 0, there is F,G e np>iL*'(0) such that for all {t, x) € [0, T] x M"^, 

\ypt{x)\ < F • (1 + \xn, \\Jt{x)\\ V \\Kt{x)\\ < G • (1 + \xfy, 

(2) there exist ei G (0, 1) and ^t G f^pyiL^i^), such that for all {t,x) G [0,r] x R'^, 

\Ao{t,x)\<^T{l + \x\^~''). 

Now we can prove 

Proposition 3.2. Under the conditions of Theorem 1.1, for almost surely w € Jig, the ODE 
(3.2) generates a unique regular Lagrangian flow Yf which leaves the Lebesgue measure quasi- 
invariant; moreover for all t > 0, Yt has a measurable inverse map Y^~ : M'^ — >■ M'^. 



Proof. We only need to check that the conditions in Theorem 2.5 are satisfied for almost 
surely w € ^q. First by (2) of Lemma 3.1, the vector field ^o satisfies the condition (1) in 
Theorem 2.5. 

For all t € [0, T], since (ft is a C^-diffeomorphism on W^ and if j : R'^ — )• W^^'^ is C^ , Lemma 
2.2 and the definition (3.1) of Aq teh us that io(i, •) G BV«oc(IR'^,IK"')- Moreover by the chain 
rule and Lemma 2.2 (1), 

DAo{t) = {DKt)Ao{^t) + KtD{Ao{^t)) 

= {DKt)Ao{^t) + \det{Jt)r^Kt [{ip^%DAo]Jt. (3.4) 

Fix any R > 0. Recall that we identify a locally integrable function / on M with the Radon 
measure f dx. Since Kt € C^, we know that {t,x) — > DKt{x) is bounded on [0,T] x B{R). 
Moreover, by the sublinear growth of Aq and Lemma 3.1 (1), it is easy to deduce the boundedness 
of AQ{ipt) on [0, T] X B{R). Hence there exists a positive constant Ct,r (depends on tu € JIq) 
such that the total variation 

\{DKt)Ao{^t)\{B{R)) < CT,RCd{B{R)), 

which implies that t —J- \{DKt)AQ{(pt)\{B{R)) G L^([0,T]). Now we consider the second term 
on the right hand side of (3.4). Again the quantities |det(Jt)|~^, Kt and Jj* are bounded on 
[0,r] X B{R), thus we only need to show that t -^ \{ip:i~^)#DAo\{B{R)) is integrable on [0,r]. 
We have 

\{^^%DAo\{B{R)) = \DAo\{y,t(.B{R))). 

By (1) of Lemma 3.1, the set UQ<t<T^t{B{R)) C B{F{\+R°^)) is bounded, from this we conclude 
that [0,r] 9 i — )■ \{lp~[ )j^DAq\{B[R)) is a bounded function, hence integrable. 

Now we check the last condition in (2) of Theorem 2.5. By Lemma 2.2 (2), we have 

|div(io(i))| < Kdiv(ifi),Ao(^i))| + |div(Ao) o(^t|. 

Since d\'v{Kt){x) and AQ{ipt{x)) are bounded on the product [0, T] x B{R), we know that t — > 
||(div(ift)5^o(<^t))||L°°{_B(ij)) is a integrable function on [0,T]. By the local boundedness of 
div(Ao) and Lemma 3.1 (1), we obtain the integrability of [0, T] 3 t ^- ||div(Ao) ° Vt\\L°°{B{R))- 
Summing up these discussions, we arrive at 



/ ||div(io(t))||L'-(iJ(i?))dt < +00. 
Jo 



Therefore all the conditions of Theorem 2.5 are verified, and we complete the proof. D 

Now we are at the position to give the 

Proof of Theorem 1.1. (Existence) We only need to show that the fiow Xt{x) = ipt{Yt{x)) 

satisfies the Stratonovich SDE (1.1). Remark that for all w a Qq ^ind any t € [0, T], the map 
Xt : M — > R is well defined almost everywhere. By the generalized Ito formula (see Theorem 
3.3.2 in [17]) and the definitions of (ft and It, we have for a.e. x € R*^, 

dXtix) = {dipt){Yt{x)) + {d^ipt){Yt{x)) dVtix) 

■ m \ 

J2 Mvt) o dwi ] {Yt{x)) + Jt{Yt{x))Kt{Yt{x))Ao [ipt{Yt{x))] dt 



:5;;^i(Xt(x)) o d^ + Ao(^t(^)) dt. 



i=\ 



To show the quasi-invariance of the flow X^, let pt be the Radon-Nikodym density of {Yt)^C(i 
with respect to Cd, then for any cj) G C^' 



cl){Xt{x))dx= l_cl)[pt{Yt{x))\dx 

<P[^t{y)Vpt{y)dy= I cp{x)[Pt\dei{Kt)\){ift\x))dx, 

where the last equality follows from the change of variable. Hence 

iXt)#Cd={pt\detiKt)\)i^^')U 

(Uniqueness) Suppose there is another solution Zf, we consider Zf = ip^ (Zt)- We will show 
that Zt solves the ODE (3.2). In fact, by [5] (see pp. 103-106, or (5.1) in [13]), 






(3.5) 



Again by the generalized Ito formula ([17] Theorem 3.3.2), 

dZt = {d^7^){Zt) + [{d^^T^){Zt)] odZt. 

Recall that for a.e. x S M , Zt{x) solves the Stratonovich SDE (1.1), and Sx-f/^^ = Ktiyf^ (a;)). 
Combining these results with (3.5) gives rise to 

Y,MZt)odw\ 



+i^i(v?,-i(Zt(x))) o Y.Ai{Zt{x))od.v^^^A^{Zt{x))dt 

Kt[Zt{x))AQ {^t[Zt{x))\ dt = io(t, Zt{x)) dt. (3.6) 



That is, Ztix) solves the ODE (3.2) for a.e. x G M . But by Proposition 3.2, this equation 
generates a unique flow If : M*^ — > M*^. Hence Yt = Zt = ^t i^t) for almost every x G M'^, which 
implies that any solution Zt to (1.1) can be expressed as the composition iptiXt)- We get the 
uniqueness of (1.1). D 

By Proposition 3.2, we know that the stochastic flow Xt in Theorem 1.1 has a inverse flow 
X^ = Yf~ o ifj which consists of measurable maps on M . 

4 Stability of (1.1) and stochastic transport equation 

In this section we study the stability of the Stratonovich SDE, proving Theorem 1.2. As an 
application we will give an explicit solution to the corresponding stochastic transport equation. 

Proposition 4.1. Suppose the conditions of Theorem 1.2. For each n>l, define AQ{t,x) := 
Kt{x)AQ[ipt[x)). Let Y^ he the unique flow generated by the ODE (3.2) with Aq being replaced 
by Aq. Then almost surely, for any T, R > 0, we have 



lim / sup |yi"(x) - Yt{x)\ dx = 0. 

^^^ Jb{R) 0<t<T 



Proof. As in Proposition 3.2, now we check the conditions in [1] Theorem 6.6. It is clear 
that condition (6.3) is verified. Remark again that the uniform boundedness assumption in 
(6.4) can be relaxed to allow uniform linear growth. Since the vector fields Aq have the same 
sublinear growth, we can prove a uniform growth estimate for AQ{t,x) similar to the one given 
in (2) of Lemma 3.1 (see [13] Lemma 2.2 for a proof). Next we check that Aq converges in 
Ll^{{0,T) X ]R'=',M'^) to io defined in (3.1). We have 

\A^{t,x)-Ao{t,x)\ < \\Kt{x)\\ . \A^{^t{x)) - Ao{M^))\. 

Since {t,x) — > Kt{x) is continuous on [0,T] x B{R), there is Ct,r > such that 

snp{\\Kt{x)\\ : {t,x) G [0,T] x B{R)} < Ct,r. 

Therefore 

[ [ \A^{t,x)-Ao{t,x)\dxdt<CT,R [ [ \A^{Mx))-Ao{Mx))\dxdt 
Jo Jb{r) Jo Jb{r) 

= Ct,r [ I 1^0 -^ol • \dei{Kt{vi^))\dxdt. 

Jo Jipt{B(R)) 

By Lemma 3.1, the set \jQ<t<T^t{J3{R)) C i?(F(l+i?")) is bounded, and the function |det(iCt(x))| 
is bounded on [0, T] x B{R). As a consequence. 



\A'^{t,x)-AQ{t,x)\dx<C'Ti{r I |^^-^o|dxdt, 

Jb(F(1+R°')) 



1 

/o Jb{r) ' Jb{f{i+r°')) 



which, by condition (2) of Theorem 1.2, converges to for almost surely w S f^o- 

Now we verify the condition (6.5) of Theorem 6.6 in [1]. By the definition of AQ{t,x), it 
holds 

VA^{t,x) = {VKt{x))A^{ipt(x)) + Ktix){VA^){^t{x))Jt{x), 

hence 

\\VA^it,x)\\ < \\\/Ktix)\\ . \A^iMx))\ + ll/^*(x)|l • \\Mx)\\ ■ ||(VA^)(^i(x))||. 

The terms ||Virt(x)||, ||i<'t(x)|| and ||Ji(x)|| are bounded on [0,T] x B{R). By Lemma 3.1 and 
the fact that Aq have uniform growth, it is easy to show that the sequence |Aq (994(3;)) | has an 
upper bound on [0, T] x B{R) independent of n. Regarding the last term, notice again that 
^o<t<T'Pt{B{R)) is a bounded subset and that |[VAq|| are locally bounded. Summing up the 
above arguments we obtain the boundedness of VAQ{t,x) on [0, T] x B{R). 

Finally to verify the condition (6.6) in [1], in view of Remark 6.3, it is sufficient to show that 
for all i? > 

n>l Jo 

By the definition of Aq and (2) of Lemma 2.2, we have 

|div(i^(t, •))! < |div(A'i)| • \A^ o ^t\ + |div(^^) o ^t\. 

Similar discussions as above lead to the desired result. Thus all the conditions in [1] Theorem 
6.6 are satisfied, and we complete the proof. D 
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Ct :=sup / ||div(i[J(t,-))|lL-{B(K))dt < +00. 

n>l JO 



Corollary 4.2. Under the conditions of Theorem 1.2, for any p > 1, almost surely 

lim I sup \Yl'{x)-Yt{x)\Pdx = fd. 

"-^°° Jb{R) 0<t<T 

Proof. By Proposition 4.1, the sequence supo<t<T |^t"'(') ~^t(')l converges to in the Lebesgue 
measure on the ball B{R). Hence we only need to show that this sequence is bounded in 
U'{B{R),dx) for any p > 1, then the desired result follows from the uniform integrability. 

By the growth estimate of Aq in Lemma 3.1, it is easy to deduce that |5^(3;)| < <I> • (1 + \x\), 
where <^ € np>iL^(rj) (see (iii) in the proof of [19] Lemma 3.2). Remark that for every n > 1, 
Aq has the same growth as Aq, hence sup„>]^ |1^"'(3;)| < <^ • (1 + |x|). Therefore 

sup |yt"(x)-yf(x)|P<2P|>P(l + |x|P), (4.1) 

0<i<T 

which implies clearly the boundedness of the sequence in L'P{B{R),dx). D 

Now we can prove Theorem 1.2. 

Proof of Theorem 1.2. We use the representations of the solutions: Xt = ftiYt) and X" = 
(pt{Y^). By the mean value formula, 

X^x) - Xt{x) = ^tiY^ix)) - ^t{Yt{x)) 

I Jt((l - u)Yt{x) + uYt^{x))dv\ {Y^ix) - Yt{x)). (4.2) 

By Lemma 3.1 and the growth estimates of Yj and 1^" given in the proof of Corollary 4.2, 

II Jt((l - u)Yt{x) + uYl'ix)) II < G(l + \Yt{x)f + \Yr{x)f) 

<G{l + 2¥)[l + \x\^). 

Therefore by (4.2), 

|Xr(x) - Xt{x)\ <(j \\Jt{{l - u)Yt{x)+uYr{x))\\d^ \Y,^{x) - Yt{x)\ 

< G(l + 2^) (1 + \x\P) \Yl'{x) - Yt{x)\. (4.3) 

As a result, for w ^VLq, 

I sup |Xj"(x)-Xi(x)|dx<G(l + 2^^)(l + |i?|^) / sup |yj"(x)-Yt(x)|dx, 

Jb{R) 0<t<T Jb{R) 0<t<T 

which converges to almost surely by Proposition 4.1. 

Now we prove the L^{il.) convergence of solutions for any p > 1. Indeed we will prove a 
stronger result: for any p > 1, 



lim E / sup |Xj"(x) - Xt{x)\Pdx = 0. 

<-^°° Jb(R) 0<t<T 



Similar to Corollary 4.2, it is enough to show that the sequence { /gf^^ supQ<j<y |Xf (x) — 
Xt{x)\P dx : n > l} is bounded in some L'^{Q,) [q > 1). However this follows easily from (4.1), 
(4.3) and the facts that G, $ belong to all L'^{rt). The proof of Theorem 1.2 is complete. D 

As an application of the stability result, now we study the corresponding stochastic transport 
equation with the purpose of constructing an explicit solution to it, by using the flow generated 
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by (1.1). First we give some remarks on the "inverse" flow associated to (1.1). To this end we 
regularize the drift vector field Aq by convolution using a standard kernel Xn'- ^o ~ 4'n'{A.o*Xn), 
here (pnix) = (p{x/n) with (p £ C^{R'^, [0, 1]) satisfying 

4>\b(i) = 1, supp((/.) C B{2). 

Let X^{x,w) be the smooth flow associated to (1.1) with Aq being replaced by Aq. Then it is 
clear that the conditions in Theorem 1.2 are satisfled, hence for any p > 1 and T,R > 0, we 
have 

lim E / sup \Xj'{x)-Xt{x)\Pdx = 0. (4.4) 

"-^°° Jb{R) 0<t<T 

Fix some T > 0, define the time-reversed Brownian motion wf = wt — wx-t and consider 

m 

dX^ = Y,M^T)°d{wfy-Ao{xT)dt, Xo = x. (4.5) 

Under the conditions of Theorem 1.1, this equation still generates a unique flow Xf , < t <T. 
Similarly we have X"' (a;, it) ) which is the solution of the above equation by replacing ^o with 
Aq, then we have (X^)~^ = X^' and 

lim E / sup \Xt''^{x)-Xj^{x)\Pdx = 0. (4.6) 

"-^°° Jb{B) 0<t<T 

With the convergence results (4.4) and (4.6) in hand, we can prove 

Lemma 4.3. Assume the conditions of Theorem. 1.1 and that the divergence div(^o) is bounded 
on W^ . Then for every T > 0, X^ = X^ a.e. on W^, and the density function ax of {X^ )#^d 
with respect to Cd has the expression: 

J2f div{Ai)iX,ix))odwi + J diviAo)iX,{x))ds 
Proof. For every n > 1, we have {{XJ^)~^):f^Cd = a^ Cd where (see Lemma 4.3.1 in [17]) 
at\x) = exp ( f; ^* div(^,)(^r(^)) ° di«: + j'^ div(^^)(X,'^(x)) d. 

Next div(^Q ) = (V(/)„, Aq * Xn) + (t>n ■ (div(^o) * Xn)- Since ^o has sublinear growth, it is obvious 
that there exists C > such that sup„>x |(^o * Xn){x)\ < C{1 + |a;|). By the definition of (j)n-, 

\{V(t>n, Aq *Xn)\< -|V0(-/n)| • |^o * Xn\ 

n 

< ^lMk(i + |x|)l|„<|,.|<2„} < 3C||V,/.|U. 

Hence the divergences 

|div(^^)| < \{V(t>n, Aq *Xn)\ + \(t>n ' (div(Ao) * Xn)\ < 3C||V,^||oo + ||div(Ao) ||oo 

are uniformly bounded on M'^. As a result (see Lemma 3.5 in [19]), for any p G R, 

sup sup sup E[(crj ^(x))^] < +00. 

n>l 0<t<T xeR'^ 

12 



Now similar arguments as in the proof of Theorem 2.6 in [23] lead to the result. D 

This lemma tells us that if div(ylo) is bounded, then the flow we constructed in Theorem 1.1 
is indeed an almost everywhere stochastic invertible flow in the sense of Definition 2.1 in [23]. 
Furthermore by (4.6), for any t > 0, it holds 

lim E /" \{Xj}y\x) - Xr\x)F dx = 0. (4.7) 

"^^"^ Jb(r] 



Now we can construct an explicit solution to the corresponding stochastic transport equation 
by using the inverse flow Xf , as in Theorem 5.1 of [13]. Though we follow the idea of the 
proof of [13] Theorem 5.1, the difference is that here we only assume that the initial value 
^0 is measurable, hence the proof of (4.19) is much more delicate than (5.18) in [13] (see [23] 
Proposition 2.3 for a different method, but 9q is supposed to be bounded there). 

Proposition 4.4. Assume the conditions in Theorem 1.1 and that div(^o) is hounded. Let 
9q :M. ^M. be a measurable function with polynomial growth. Then 9{t,x) := Oo{X^ (x)) is a 
distributional solution to the following stochastic transport equation: 

m 

d9{t) = -Y,{^m,A^)odwi-{V9{t),Ao)dt, 9\t=o = 9o, (4.8) 

that is, for any (p G C;?°(M^), 

m rt rt 

i9{t), c^)l2 = (90, 4>)l2 + V / (0(s), div(</.A,))L2 odwl+ [9{s), div(c/.^o))L2ds. (4.9) 

i=i Jo Jo 

where (•, •)^2 is the inner product in L^(]R , dx). 

Proof. First we transform the equation (4.9) into the Ito form: 

(0(t),0)^2 = (0o,</')L2+V / (^(s),div((/)A,))i2du;:+ / {9{s)AH^Ao))L^ds 

i=i Jo Jo 

+ 2^/ (^(s),div(div(</)^i)^i))L2ds. (4.10) 

The proof is similar to that of Theorem 5.1 in [13] and we divide it into two parts. 

Step 1. We assume ^o £ C"^ . In this step, similar to the proof of Theorem 5.1 in [13], the 
key point is to show 



lim E sup \9n{t,x)-9{t,x)\'P\ =0, 
"^~ \\x\<R J 

where 9n{t,x) = 9o[{XJ^)~^{x)^ and X" is defined before Lemma 4.3. However, similar to (5.12) 
on p. 1102 of [13], we only need a weaker form of the convergence result: 

lim E / \9nit,x)-9{t,x)\fdx = 0. (4.11) 

"■^"^ Jb(R) 



To this end, notice that \9n{t,x) - 9{t,x)\ < \\V9o\\oo\{XJ})~^{x) - X;''^{x)\, thus by (4.7), we 
still have 

E / \9n{t,x)-9{t,x)\Pdx< llV^ollSoE / \{X]})-\x) - X^^{x)\Pdx ^ 0, 
Jb{r) Jb{b) 
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as n — 7- oo. Hence (4.11) holds and the rest of the arguments in Step 1 of the proof of Theorem 
5.1 in [13] stih work here. 

Step 2. Now suppose that ^o is measurable with polynomial growth. Define 6q = (pn (^o*Xn)- 
Then there exist C > and qo G Z+ which are independent of n such that 

\Oo{x)\ < C (1 + |x|'?o), X G M'^. (4.12) 

Use again the notation On{t, x) to denote 6n{t, x) = 6q{X^ {x)) where Xt{x) is now the solution 
to SDE (1.1). Then by Step 1, On satisfies (4.10). Now using the SDE (4.5) and the moment 
estimate, we have for any T > and t G [0, T], 

e( sup \Xl{x)\A <Cp^T{l + \x\^)- 

In particular, for each t G [0,r], 

E(|xr^(x)n<c,,T(i + ixn. 

By (4.12), it holds that 

supE(|6'„(i,x)|P)+supE(|0(t,2;)|P) < Cp,r (1 + I^P). (4.13) 

t<T t<T 



Therefore for any p>2 and (/> G Co°(M'^) with supp((/)) C B{R) 



[ E{\{enit)A)L^\'')dt<( [ \4>\'^dx) II E(|0„(t,x)r)dxdi 

JO \JK.d J Jo Jb(r) 

<CpT\ I Wdx] I (l + |xP)dx<+cx). 
' Vjm'' / Jb(r) 



(4.14) 



where q is the conjugate number of p. 
Fix some M > 0. We have 



e/ \en{t,x)-e{t,x)\^dx = W. I \e^{y)-eo{y)\^at{y)dy 

Jb{r) Jx-\b{b)) 

<^( I +1 )K{y)-Oo{y)\^at{y)dy 

\Jb(M) J X-^ (B(R))\B(M) J 



B(M) JX-\B{R))\B[M) 

im+m)- (4.15) 



Since div(Ao) is bounded, we have by Lemma 3.5 in [19] that 

sup sup E((Jt(x)) < Ct < +00, 

0<t<Txi 



hence 

limsup/r(t) =limsup / \9^{y) - 9o{y)Mat{y)) dy 

n->oo n— i-oo J B(M) 

<Ct lim I \e^iy)-eo{y)\^dy = 0, (4.16) 

"-^°° JBiM) 

where the last equality follows from the convergence of 9^ to 9q in Lf^^{'R.'^). By the Cauchy 
inequality. 



I^{t) <U I K{y) - eo{y)\%{y)dy) (e / ^t{y)dy) 

V Jx-\b(r))\b(m) J V Jx-\b(r))\b(m) J 
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=--{l^,iit)Ih(t)f^'- (4-17) 

We have by (4.13), 

/£i(t)<E/ K{y)-Ooiy)\%{y)dy 

JX-\B{R)) 

= E [ \9^iXr\x))-9o{Xf\x))\'dx 
Jb{r) 

<C' I (l + |3;|^««)dx<+oo. (4.18) 

Jb{r) 

Next the function ^t'^x~^(B{R])\B(M) ^^nds to as M tends to +oo for P x £^-a.e. {w,y) € 
fio X M'^; moreover (yax-\B{R))\B{M) ^ ^t'^x-\B[R)) ^nd 

E / ^ atlx-\B[R)) ^y = ^ j ^ lB(R)dy = Cd{B{R)) < +(^. 
Hence by the dominated convergence theorem, we have 

This plus (4.17) and (4.18) tells us that 

lim mt) = 0. 

Therefore by (4.16), first letting n goes to +oo in (4.15), and then letting M goes to oo, we 
obtain 



lim E \9n{t,x)-e{t,x)\^dx = 0. 

'^+~ Jb(R) 



From this we deduce that for any (p £ C^{W^) with supp(<^) C B{R), 

E[{{9n{t),<P)L2-{e{t),4>)L2)^]<( f <P^dx)E [ \en{t,x) -9{t, x)fdx^O (4.19) 

V Jr'' j jb(r) 

as n ^^ CO. Now (4.14) and (4.19) allow us to pass to the limit and the proof is complete. D 

Now we discuss the connection between the stochastic transport equation (4.8) and the 
following transport equation associated to the random vector field ^o defined in (3.1): 

dut = -{Vut,Ao{t))dt, u\t=o = uo. (4.20) 

To this end we first give some preparations. Recall that (ft is the smooth flow defined at 

the beginning of Section 3. Let pt = '^^^'^^l*^"^ and pt = '^^^%'c*^''^ ^e the Radon-Nikodym 
densities. We have the following simple equality: 

pt{x) = [p,{^^\x))]-\ (4.21) 

Indeed, for any ip G C^(M'^), we have 

i^[^t{y)]pt{y) dy = / Tp{x)pt{ip^^{x))pt{x)dx, 
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which leads to (4.21) due to the arbitrariness oi ip £ C^{M.'^). Moreover by Lemma 4.3.1 in [17], 
the density pt has an exphcit expression: 



pt{x) =exp(f2 [ dw{Ai){ips{x)) o dw'^ . (4.22) 



Now we show that the distributional solutions of (4.8) and (4.20) are related to each other 
by the smooth flow ipf 

Proposition 4.5. Suppose that 9t is a distributional solution to the stochastic transport equation 
(4.8), then almost surely, Ut := Ot{'^t) solves the transport equation (4.20) with u\t=Q = 9q. 

Proof. For any f/' G C^(M'^), we have 

0i(v?i)V'dx= / eti^{^^^)ptdx. (4.23) 



Now we compute the Stratonovich stochastic differentials of V'(¥'t ) ^^"^ Pt- ^Y [5] (see pp. 
103-106, or (5.1) in [13]), 

m 

dc^^i = -Kt{^^^)Y,Mx)odwl. (4.24) 

Hence 

m 

d^((/.,-') = ((VV')((^r'),°d(^r'> = -Y,{{Vi'){^i^),Kt{^i^)A,)odw\. 

i=l 

Notice that V{tp{Lp~[^)) = Kl{ip^^){Vip){ip^^), we obtain 



dV'(^r') = - E (V(V'(^r')), A^) ° M- (4.25) 

Next we compute dpt- By (4.22), 

m 

dpt = Pt E 'iiv(Ai )((/?() o dw\, 

i=l 

hence we deduce from (4.24) and the generalized Ito formula that 

d[pt{^i^)\={dpt){^i^) + {{Vpt){^i'),od^-^) 
m m 

=pt(^r')Edi^(^*) °d^* -^{{K:vpt){v^'),A,)odwi. 

i=l i=l 

Using again the Ito formula and by the relation (4.21), we arrive at 
dpt = - [pti^fT^)] ^^ o d[pt{ip^'^)] 

m m 

= -Pt E di^(^^) ° d"^t + Pt E {^iptivt^))^^i) ° M- 

i=l 1=1 

Since Vpt = —plV{pt{iPi )), finally we obtain 

mm m 

d/Ot = -Pt E div(^i) o dwi - E(VPt' ^i) ° d«ij = - E div(A^i) ° du;j. (4.26) 

4=1 i = l 4 = 1 
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Now by the equalities (4.25), (4.26) and the fact that 6^ solves the stochastic transport 
equation (4.8), we have 



d[dti;{^^')pt]=^{ipT')pti-Y,{V9t,A,)odwi-{^et,Ao)dt 

m m 

i=l i=l 

Since (li\{4i{ip~[^)ptAi) = ip{ip'^^)(ii\{ptAi) + pt{V{ip{(pt'^)),Ai), we arrive at 



d[eti;{^i^)pt] = HipT')pt { - Y,{V9t,A,) o dwi - {VOu Ao) dt 

m 

- ^ 6itdiv(^(v?7^)ptAi) o dwi. 



i=l 



The above equality should be understood in the distributional sense. More precisely we have 
obtained 



/ etil^{ip^^)ptdx= ^oV^dx + V/ ( / e,d\Y{'4){ip,^)psA,)dx]oduf, 

JR'' jR'i ^^^ Jo \ JR'i J 

+ / / 0sdiv(V'(v37^)p,^o)dxds 

JoJR<^ 

-Jlj (j^Jsdiv{^l;iip-')psA,)dx]odwi 
= / 6'oV'dx+ / / esdiv{4^{ip-'^)psAo)dxds. 

jRd JoJRd- 

As in Proposition 4.4 we denote by (•, •)^2 the inner product in L"^ {W^ , dx) . By (4.23) we have 

{et{ipt),iP)L2 = iOo,iP)L2 + [ I esdiv{i;{ip;^)psAo)dxds. (4.27) 

JoJr'^ 

By the definition of ps, 

sdiv{ip{ipJ^)psAo)dx 
= [ 9sPs[{iKyi;)iipj'),Ao) +i;i^;^)p;\Vps,Ao) +Hv;')div{Ao)]dx 

jRd- 

= / esiips)[{'^i^, KsAoiifs)) + ij{ipJ^V ps)iips), Ao{ips)) + il^div{Ao){ip,)]dx. 
Jr'' 

Lemma 2.3 leads to 

div(io(s)) = (div(ir,),Ao(¥?s)) +div(^o)('Ps) 

= ((P7^Vp,)((^,), Ao{ip,)) + dw{Ao){ips), 

it follows that 

6sdw['il^{ip~^)psAo)dx = / 9s{(ps)div{'ipAo{s))dx. 



i2ds. 



This plus (4.27) gives us 



Jo 
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which means that almost surely, Of{ipt) is a distributional solution to the transport equation 
(4.20) with initial value ^o- Q 

Remark 4.6. Originally we intended to prove the uniqueness of the solutions to the stochastic 
transport equation (4.8) by using the above proposition. Indeed, if any solution of (4.20) can 
be represented as ut = ^0(1^^" ), where Yj is the flow generated by Aq, then by the above 
proposition, we must have Otift) = ^o(^i~ ) (since ^t((/'t)|t=o = ^o); which gives us 

That is to say, any solution of (4.8) is expressed as the composition of ^0 and the inverse flow 
Xf . However, since the divergence div(^o) of ^0 is unbounded, it is difficult to get a meaningful 
uniqueness result for the equation (4.20), see [1, 2, 9]. 

5 Approximate differentiability of the flow generated by (11) 

In this section we study the approximate differentiability of the stochastic flow X^ associated to 
the Stratonovich SDE (1.1) whose drift coefficient Aq belongs to the Sobolev space Wi^^ (R'^,M.'^). 
To this end, we introduce some notations and results about maximal functions. For any bounded 
measurable subset U C M. with positive Lebesgue measure Cd{U) > 0, define the average of 
/ G LJjR'i) on U by 

-/ f{x)dx = ——— / f{x)dx. 
Then for any x G M*^ and R > 0, the local maximal function Mjif is defined by 



MRf{x)= sup / |/(y)|dy 

0<r<i? JB(x,r) 



{x,r) 

where B{x,r) = {y (zW^ : \y — x\ < r}. Here are some results regarding the maximal function 
(see [22]; for a proof of the second result, cf. the Appendix of [15]). 

Lemma 5.1. (1) For R, p > 0, there are C^, C^^p > such that for all f G -L;^^^(R"'), we have 
[ Mnf{x)dx<Cd,p + Cd [ |/(x)|log(2 + |/(x)|)dx 

JB{p) Jb(R+p) 

and for any a > 0, 



Caix G B{p) : Mnfix) >a)<^ [ |/(x)| dx. 

« Jb(R+p) 



B{R+p) 

(2) Let f G W^^^ (M'^) . Then there is Cd > (independent of f ) and a negligible set N C W^, 
such that for all x, y G N'^ with \x — y\ < R, 

|/(x) - fiy)\ < Cd\x - y\{{MR\Vf\)ix) + (MR|V/|)(y)). 

We first prove the following result on the approximate differentiability of the regular La- 
grangian fiow generated by a Sobolev vector field b. This is an extension of Corollary 2.5 in [7] 
to the case where b has linear growth (see [7] Corollary 3.5 and [4] Remark 3.8 for more general 
case, but therein the divergence of b is assumed to be bounded on M ). 

Proposition 5.2. Assume that b G L^{[0,T],wl;^{R'^,R'^)) satisfying 
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I''t(^)l c roo/rn t^i ^ md 



(i) m£jeLooQo^j.] 



X 



pT 



(ii) for any R>0, f^^ ||div(6j)||ioo(B(R))dt < +00; 



(iii) for any R>0, J^J^iR) 1^^*! ^°s(2 + |V6t|) dxdt < +cx). 

Lei Yi be the regular Lagrangian flow generated by b. Then for any R > and e > 0, there exists 
a Borel set E C B{R) such that Cd,{B{R) \E) < e and the restriction Yt\E is a Lipschitz map 
for any t G [0,T]. 

In particular, for any t € [0, T], If is approximately differentiable Cd-a.e. in W^ . 

Proof. We follow the idea of the proof of Corollary 2.5 in [7]. For 0<t<r, 0<r<2i? and 
xe B{R), define 

Q{t, x,r) = f log h 1 dy. 

JB{x,r) \ r J 

Prom Definition 2.4(1), it follows that for a.e. x and for all r € (0,2i?], the map t — > Q{t,x,r) 
is Lipschitz and 



-T-(i,x,r)< -/ 
dt Jb(xa 



dt ^ ^ dt ^^' 



■{\Yt{x)-Yt{y)\+r) Vy 
bt{Yt{x))-bt{Yt{y))\ 



dy. (5.1) 



lB(x,r) \Yt{x)-Yt{y)\+r 
By condition (i) and Gronwall's inequality, it is easy to show that 

\Yt{x)\ < (1 + R)e'^'^, for ah x e B{R), 0<t<T. 
Therefore for a.e. x G B{R) and y G B{x, r), we have 

\Ytix) - Ytiy)\ < \Ytix)\ + \Ytiy)\ < (1 + R)e^^ + (1 + 3R)e^^ = 2(1 + 2R)e^^ =: R. 
Since {Yt)#^d ^ ^d-, we can apply Lemma 5.1(2) to get 

\bt{Yt{x)) - bt{Yt{y))\ < Cd\Yt{x) - Yt{y)\ • [{M^\Vbt\){Yt{x)) + (M^|V6f |)(y,(y))] . 
Substituting this estimate into (5.1) gives us 

^{t,x,r) < / Cd[{M^\Vbt\){Yt{x)) + {M^\Vbt\){Yt{y))]dy 

a* JB(x,r) 

= Cd{Mj,\Vbt\){Yt{x)) + Cd / {Mf,\Vbt\){Yt{y))dy. 

JB(x,r) 



J(x,r) 

Integrating with respect to time, we obtain for all t € [0,T], 

cT rT 

?(a;,r) 



Q{t,x,r) <\og2 + Cd I {M^\Vbs\){Ys{x))ds + Cd f / {M^\Vbs\){Ys{y))dyd, 

Jo Jo JB(x,r) 



Let $(x) = J^ {Mp^\Vbs\){Ys{x)) ds, then by Fubini's theorem. 



Q(i,x,r) < log 2 + C7rf$(x) + Crf / Hy)dy, for alH G [0, T]. 

JB{x,r) 
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Hence by the definition of the maximal function, 

sup sup Q{t,x,r)<log2 + Cd^{x) + Cd{M2R^){x). (5.2) 

0<t<T 0<r<2R 

For 1] sufficiently small, we have 

Cd{x G B{R) : log2 + Cd<^{x) + Crf(M2R$)(x) > 1/r/) 

<Cd{xe B{R) : CdHx) > l/(3r?)) + Cd{x £ B{R) : Cd{M2R^){x) > 1/(37/)) . (5.3) 

By Chebyshev's inequality, 

Cd{x G B{R) : Crf$(x) > 1/(37?)) < 3i]Cd [ $(x) dx. 

Jb{r) 

Using Lemma 5.1(1), we have 

Cd{x € B{R) : Cd{M2R^)ix) > 1/(37?)) < 3r]CdC'd [ Hx) dx. 

Jb{3R) 

Substituting these two estimates into (5.3) and by the definition of ^{x), we obtain 

/ := Cd{x € B{R) : log 2 + CdHx) + Cd{M2R<^){x) > I/7?) 

< 3r]Cd{l + C'd) [ Hx) dx 
Jb{3B) 

= SrjCdil + C'a) [ [ {Mj^\Vbt\){Yt{x)) dxdt. 

Jo Jb(3R) 



'0 JB(3R) 

Using the density pt of the flow Yj, we get 

/ < 3vCd{l + C'a) [ [ {Mj^\Vbt\){y)pt{y) dydt. 

Jo JYt{B{3R)) 

In view of the expression of pt given in Remark 2.6, for any x € B{3R) and t G [0, T], 

Pt(Yt{x)) = exp f - / div{bs){Ys{x))dsj < exp f / ||div(6^)||i<x=(B(i?i))ds j =: L, 

where Ri = (1 + 3R)e . Hence by Lemma 5.1(1), 

/ < 3vCd{l + C'a)L f f (M^|V6i|)(y) dydt 
Jo Jb{Ri) 



< 3rjCd{l + C'a)L [ Cd,R, + C'^ f |V6t| log(2 + |V6t|) dy 

Jo L Jb{Ri+r) 

=:7?Li. 

Now for any e > 0, let r? = e/Li, then by (5.2) and the definition of /, we have 
Cdi X & B{R) : sup sup Q{t,x,r) > — ) < I < ^ ■ Li = e. 

V 0<t<T 0<r<2R £ / ^1 

Let 

E = <x £ B{R) : sup sup Q{t, x,r) < — > , 

I 0<t<T0<r<2_R e J 
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dt 



then j:d.{B{R)\E) < e and ior any x e E, < t < T and < r < 2R, the definition of Q{t, x, r) 
leads to 

u\xM^m\^,U<^^ (5.4) 

B{x,r) \ r J £ 

Now fix any x,y (z E and let r = \x — y\ which is less than 2R. We have by the triangular 
inequality, 



JBi 



i„^ , y.W - YM\ ^ \ ^ ^^ / |y,(x) - y,(.)| + \Y,(z) - Yfy)] ^ ^ 



<iog|l^'W-^''--)l+iViog('l^''--)-^-fa'l+i 



therefore by (5.4), 

JB{x,r)nB{y,r) \ r J 

?(x,r)nB(j/,r) 






JB{x,t) \ r J JB{y,r) 

e 
where C^ = Cci{B{x, r))/ Cii{B{x, r) n B{y, r)) only depends on the dimension d. Therefore 

\Yt{x) - Yt{y)\ < re^^^^^l^ = \x- y\(?^^^^l^ 

which holds for all x,y e E. Hence Lip(Yt\E) < e^^'^^^l'^ . D 

Now we can prove the main result of this section. 
Theorem 5.3. Assume that Ai, • • • , An € C^+'^(M'^,M'*) and that Aq G Wl;^{R'^,R'^) satisfies 

(1) Aq has suhlinear growth; 

(2) div(ylo) is locally bounded on W^; 

(3) for any R> 0, J^^j^^ ||VAo|| log(2 + ||V^oll) dx < +oo. 

Then for a.s. w € ^o, for any R > and 5 > 0, there exists a Borel set E C B{R) such that 
Cd{B{R) \E) < 6 and the restriction of the flow Xt to E is a Lipschitz map for any t € [0,T]. 
In particular, Xt is approximately differentiahle Ld-a.e. in W^ for any t G [0,T]. 

Proof. Since Xt = V'f(^) ^-^id for a.s. w € i^O) the map (/?( : M^ ^ M'^ is a C^-diffeomorphism 
on M , we only have to prove the result for the solution Yt of the random ODE (3.2). Now we 
check that Aq satisfies the conditions given in Proposition 5.2. First by the definition of ^o(^) ■) 
and the conditions on Aq, it is clear that io(i, •) e E^{\^,T\,L]^J^^ ,W^)). Moreover 

Vio(i,x) = (VKi(a;))Ao(c/^i(x)) + A-i(x)(VAo)((/^i(x))Ji(x), 

hence 

||Vio(t,^)|| < ||Vi^t(x)|| ■ \AQ[^t[x))\ + ||i^i(x)|| • \\Jt{x)\\ ■ \\{VAo){^t{xm. (5.5) 
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The terms ||Vi^i(x)||, ||i<'i(x)|| and ||J((a;)|| are bounded on [0,r] x B[R). By Lemma 3.1 and 
the fact that Aq has subhnear growth, it is easy to show that |Ao((/?f(x))| has an upper bound 
on [0,r] X B{R). As for the last term in (5.5), noticing that L := DQ<t<T'^tiB{R)) is a bounded 
subset, we have 



// \\VAo{t,x)\\dxdt<CT,R + C^n [ [ \\{VAo){ipt)\\dxdt 

Jo J B{R) ' Jo Jb{R) 

= Ct,r + C'tr I I llVAoll • \det{Kt){^i^)\dxdt 

10 Jipt{B{R)) 



< Ct,r + Ct^r / llVAoll dx < +00, 



where the last inequality follows from the boundedness of |det(ifj)| on [0,T] x B{R). Hence 

AoeL^{[o,T],wl;liM.'',R'')). 

By Lemma 3.1, the condition (i) in Proposition 5.2 is easily checked for ^o- The second 
condition (ii) can be verified by using the equality in Lemma 2.2(2), as we have done at the end 
of the proof of Proposition 3.2. 

Now we check that ^o satisfies the condition in Proposition 5.2(iii). Again by (5.5) and 
the above discussions, we have ||VAo(t)|| < Ct,r{^ + ||(V^o) ° '/'tlD- Therefore by the simple 
inequality log(l + s) < s for all s > 0, we have 

log (2 + llVio(i)ll) < log [(2 + Ct,r){2 + \\{VAo) o ^t\\)] 
< (1 + Ct,r) + log (2 + ll(VAo) o 9^* ID- 
AS a result, 

||Vio(i)||log(2+||Vio(t)||) 

< Ct,r(1 + Ct,r){1 + ll(VAo) o ^t\\) [I + log (2 + ll(V^o) o ^t\\)] 

< Ct,/?(1 + Ct,r) [2(1 + ll(V^o) o iptW) + ll(VAo) o ip,\\ log (2 + \\{VAo) o (^J)]. (5.6) 

Again by the fact that L := Uo<t<TftiB{R)) is bounded for any R > 0,we have by the condition 
(3) that 



/ / ||(V^o)o'/'t||log(2 + ||(V.4o)o^i||)dxdi 
Jo Jb(r) 



< [ [ ||V^o||log(2+||V^o||)|det(/Ct)o^-i|dydt 

Jo Jipt{B{R)) 

< C't^rT j \\VA4 log (2 + ||V^o||)]dy < +oo, 

where C^^^ = sup{|det(is:f(x))| : (i,x) G [0,T] x B{R)} < +oo. This and (5.6) clearly imply 
that ^ 

n||Vio(i,a;)||log (2 + \\SJAQ{t,x)\\)dxdt < +cx3. 
_ J{R) 

The condition (iii) in Proposition 5.2 is verified and the proof is complete. D 
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